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Abstract 

Using the Kadanoff-Baym non-equilibrium Green's function formalism, we 
derive the self-consistent Hartree-Fock-Bogoliubov (HFB) collisionless kinetic 
equations and the associated equation of motion for the condensate wavefunc- 
tion for a trapped Bose-condensed gas. Our work generalizes earlier work by 
Kane and Kadanoff (KK) for a uniform Bose gas. We include the off-diagonal 
(anomalous) pair correlations, and thus we have to introduce an off-diagonal 
distribution function in addition to the normal (diagonal) distribution func- 
tion. This results in two coupled kinetic equations. If the off-diagonal distri- 
bution function can be neglected as a higher-order contribution, we obtain the 
semi-classical kinetic equation recently used by Zaremba, Griffin and Nikuni 
(based on the simpler Popov approximation). We discuss the static local 
equilibrium solution of our coupled HFB kinetic equations within the semi- 
classical approximation. We also verify that a solution is the rigid in-phase 
oscillation of the equilibrium condensate and non-condensate density profiles, 
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oscillating with the trap frequency. 
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I. INTRODUCTION 



In their classic book, Kadanoff and Baym (KB) ||] developed a systematic way of de- 
riving quantum Boltzmann equations using non-equilibrium Green's functions. KB derived 
generalized Boltzmann equations for interacting Bose and Fermi gases in the normal phase 
(no broken symmetry). The advantage of this approach is that one can generate kinetic 
equations starting from a well-defined single-particle self-energies using functional deriva- 
tives. The KB method was generalized by Kane and Kadanoff (KK) to deal with a 
Bose-condensed gas, with the goal of using the resulting kinetic equations to derive the 
two-fluid hydrodynamic equations of Landau. The KK analysis worked with diagonal and 
off-diagonal self-energies, and included terms of second order in the two-particle interac- 
tion (the Beliaev second order self-energy diagrams 0). This gives kinetic equations which 
include the effect of collisions between the atoms. 

In the present paper, we use the KK approach to deal with a trapped Bose-condensed 
gas. For simplicity, in the present paper we limit ourselves to Hartree-Fock (HF) self-energies 
and thus the kinetic equations we obtain are in the collisionless approximation Q. Because 
we include both normal and anomalous pair correlations, our self-energies give the so-called 
Hartree-Fock-Bogoliubov (HFB) approximation, as reviewed by Griffin |Q. While the HFB 
approximation in a uniform gas is known to lead to a single-particle spectrum with an energy 
gap in the long wavelength limit, it is a "conserving" approximation which will generate two- 
particle response functions which satisfy conservation laws (see Refs. |l]J^] and Section VI 
for further discussion). 

Our main formal result is given in Section III, namely two coupled kinetic equations 
for the diagonal (fi) and off-diagonal (/ 2 ) distribution functions for the excited atoms, in 
addition to an equation of motion for the condensate order parameter. If f 2 can be ignored, 
we obtain the semi-classical collisionless kinetic equation for fi which has been used in recent 
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discussions of the non-equilibrium properties of the non- condensed atoms |],[7) . In order to 
gain some insight into our HFB kinetic equations, in Section IV we discuss the static local 
equilibrium solutions within the semi-classical approximation and verify that these satisfy 
our equations. In section V, we show that our equations have a solution which corresponds 
to a a rigid oscillation of the equilibrium density profiles (normal and anomalous), with a 
frequency equal to the parabolic trap frequency This generalized Kohn mode is generic but 
how it arises as a solution of our HFB kinetic equations gives one insight into their structure. 



II. DERIVATION OF HFB EQUATIONS 

In terms of Bose quantum field operators, the many-body Hamiltonian (K = H — fioN) 
describing interacting Bosons confined by an external potential U ex t(r) is given by: 



K 



drift' (r) 



2m 



i[)(r) + — / drdv' ijj\v)ijj\v')v (r — r')ip(r)ip(r'). 



(1) 



For a discussion of the properties of a dilute Bose gas at low temperatures the two-body 
interaction v (r — r') can be effectively treated using the s-wave approximation: v(r) = gS(r), 
where g = Aira/m (we set h = 1 throughout this article). We shall be interested in the 2x2 
matrix single-particle Green's function defined by p]I| 



g(l } l'-U) = -i 



1 (T^(l)^t(i')) (1Y(1)V(10> ^ 



V 



(2) 



Here, T represents the time-ordering operation and we use the usual abbreviated notation 
1= (rt) and V = (r't'). We separate out the condensate part of the field operator in the 
usual fashion 



#0 = (V(r)) t + V(r), 



(3) 
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where (0( r )) = an d (V'( r ))t = ^( r ;0 is the macroscopic wavefunction. The non- 
condensate (or excited atom component) field operators ip(r) and if)^(r) satisfy the usual 
Bose commutation relations. Using (3), the matrix in (2) splits into two parts 



g{l,l';U)=${l,l';U) + h(l,l';U). 



(4) 



Here g is identical to (2) except that it involves the non-condensate part of the field operators, 
and 



MM'; U) 



( 



'mm \ 



$(i)$*(i') $(i)$(i') 



(5) 



with (ifj^(r)) t = $*(r,£). In a Bose-condensed system, the finite value of $(r,t) leads to 
finite values of the off-diagonal (or anomalous) propagators and (l)^^ (!'))• 

These must be dealt with on a equal basis with the diagonal (or normal) propagators, which 
is the reason we must work with a 2 x 2 matrix single-particle Green's function. 

A very convenient and elegant way of generating the equations of motion for g and h is 
to use functional derivatives with respect to weak external fields. The latter are described 
by 

H'(t) = J dr[U{r, t)n(r, t) + r,(r, t)ft(r, t) + r)*(r, t)$(r, t)}, (6) 

where U(r,t) is the external generating scalar field, while t](r,t) and r]*(r,t) describe the 
symmetry-breaking fields involving "particle sources" PJlQfl. The higher-order Green's func- 
tions can all be expressed as a functional derivatives of single-particle Green's functions with 
respect to these fields. 

Following the Kane-Kadanoff (KK) analysis fjJH, the HFB equations of motion can be 
conveniently written in the following 2x2 matrix form 



dl 



$o _1 (U) 



(1,1) |(1,1') = 5(1,1'), 



(7) 



^(l, 1') - t HF (l, 1') =5(1,1') 



(8) 



dl 



^(l, 1) " ^(1, 1)1 1') = -^(^(^(l')), 



dlh(l,l) g^(l, l')-^(M') = -i(*(l)>^(l'). 



(9) 



(10) 



In the above equations, integration over dl means integration over the coordinates (rj,ii), 
5(1, 1') = 5(r — r')S(t — t'), and we have to introduce spinors 

( ^ 



(*(!)> 



, (^t(i')) = (*•(!), $(!)) 



V $*(l) y 



(11) 



Here r} ea;t describes the external particle-source fields defined in (6), 

f) ext {l) 



( 77(1) ^ 



(12) 



\v m Q) j 

The inverse of the 2x2 matrix non-interacting propagator go{l, 1') is defined by 



IT 



{3) — + P^-U ext (T 1 )-U(l)+H> 



5(1,1')- 



dti 2m 

In the HFB approximation, the 2x2 self-energies in (7)-(10) are given by 

1 2n(l), m(l) \ ( 2n(l) +njl), m(l) \ 

v m*(l), 2n(l), 



(13) 



S HF (1,1')=^ 



5(1,1'), S HF (lA')=g 



' 2ra(l) + ra c (l), mil) 
v m*(l), 2n(l)+n c (l) j 



5(1, 1'). 
(14) 



In the above equations, n(l) is the total local density given by 



n{\) = i \~g<(l, 1+) + fc u (l, 1)1 = $t(i)^(i))+ | $(l) |2= ^(i) + 7lc (i) > 



(15) 



where n(l) and n c (l) are non-condensate and condensate density, respectively Similarly, 
m(l) is the "anomalous" local density defined by 
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m(l) = i [~g<(l, 1) + h 12 (l, 1)] = (0(1)0(1)) + Ml)} 2 = m(l) + [$(l)f 

m*(l) = • 1) + Ml, !)] = (^(1)^(1)) + r (I)] 2 = m*(l) + [<F(l)f . (16) 

These HFB results were first written down in this formalism by Kane and Kadanoff 0, 
with the "particle-source" fields rj and rj* left implicit. The equation of motion for the order 
parameter $(1) in this HFB approximation is given by 



1 J" + P~ + ^ ~ C/ -*( r i) " U W ~ 9 [2n(l) + n c (l)] 



$(l)=^m(l)$*(l)+r ? (l), (17) 



and its complex conjugate. Equations (7)-(0) are a closed set of equations and define what 
is called the dynamic HFB approximation. 

We now turn to solving the HFB equations of motion for atoms in the presence of a 
trapping potential. If the external generating fields induce a disturbance with a wavelength 
much longer than thermal wavelengths and frequencies much smaller than characteristic 
particle energies then the propagator g(l, 1') = g(r,t; R, T) can be expected to vary slowly 
as a function of the center-of-mass coordinates 

R=^(r x + !•!,), T=±(t 1 +t v ), (18) 

and to be dominated by small values of the relative coordinates 

r = n-n», t = ^[t x -t v ). (19) 

More precisely, Fourier transforming with respect to r and t, the function g(p,uj;H,T) 
describes the density of elementary excitations of momenta p and energy u at point (R, T) 
Jl|. These quasiparticles are assumed to have high momentum and energy (relative to the 
collective modes we would obtain from the kinetic equations), which means that g(r, t; R, T) 
is mainly weighted at small values of the coordinates r and t. 

If the Bose-condensate order parameter (0(1)) is written in terms of amplitude and phase 
variables 



(20) 



we can generalize the usual definitions for superfluid velocity and local chemical potential 
to non-equilibrium systems by identifying p| 



Vi0(l) 
89(1) 
dt 1 
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-fi + -mv s 2 {l) 



(21) 



The superfluid velocity v s enters as the gradient of the phase of the condensate wavefunction 
and the local chemical potential /x(l) is connected with the time-derivative of the phase. 
However, a problem arises in that the phase is a rapidly varying function of (R, T) which 
induces strong variations in the off-diagonal elements of h, and these are coupled to the 
components of g. To remove this strong (R, T)-dependence associated with the phase 9, we 
apply the well-known || gauge transformations on h(l, V) and g(l, 1'): 



h>(l,l') = e- ie ^ i3) h(l,l')e i9 W Ti3) 



(3) 



(22) 



where is the Pauli spin matrix. The physical interpretation of ( P^D is that it involves 
a transformation to a coordinate system in which non-condensate atoms are moving with 
average velocity v s with respect to a stationary condensate. The gauge transformation ( |22|) 
removes the strong (R, T)-dependence associated with the order parameter phase 9 and 
leaves the equations of motion (7)-(10) invariant if we replace g^ 1 in (13) by 



9o 



IT 



(3) 



5(1,1')- 

(23) 



After carrying out this gauge transformation, the HFB equations for gf^l, 1'; U) and 

< 

<7f^(l, 1'; U) given by (7)-(10) become 
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09(1) 






% 0t x 


Of, 


^ 2m 

— III 




' . 


09(1) 








dt 1 


+ 2m' 


. 


99(1') 


-iv 


% 0t v 


dt v 


". 


09(1') 

Oty 


2m l 


l dt v 



[Vl+zV^(l)] 2 + /io-f/e//(l^ 



< 
~,> 



< 
~,> 



i'M2 



i 'M2 



^(M';[0=W(%(M'; 

< < 

</f 2 (l,l';£7) = ^'(1)^22(1, 1'; £7) 
</>(l, l';£7) = ^(10^12(1, i';E0 

</>(!, 1'; C7)= ^'(l'V^l, 1'; CO- (24) 



Here 



t/ e// (l) = U ext {r x ) + £7(1) + 2gn'(l) (25) 
is the effective self-consistent Hartree-Fock dynamic mean field. The condensate part 2<m c (l) 



in ( 25]) can be viewed as an additional "external field" acting on the non-condensate. Since 
we will work with these gauge-transformed correlation functions in the rest of this paper, 
we drop the primes on gu, g\2, n and m to simplify the notation. 



The corresponding equation of motion for the condensate amplitude Jn c (l) is in the 
moving frame of reference can be written in the form: 



° ^ 1 ] Vl \mw 2 s + /i - U ext (n) - £7(1) - g [2n(l) + n c (l)\ + 



Oh Oh 



2m 



+iv s (i)- Vi + -Vi-v s (i; 



V^ = ^m(l) v Ql) + r/'(l), 



(26) 



where T]'(l) = rj{l)e~ ie ^ 1 ' is the symmetry-breaking source function in the moving frame of 
reference. We note that there is a factor of two difference between how the condensate and 
non-condensate enter these equations. This is because atoms in the condensate are in the 
same state and thus there is no exchange part. In the case of non-condensate atoms, both 
Hartree and Fock terms arise since we are dealing with atoms in different states. 
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III. HFB KINETIC EQUATIONS 



To rewrite the equations of motion derived in Section II in the form of kinetic equations, 
we recall the connection between the usual single-particle distribution function /i(p, R, T) 
and the diagonal Green's function g^ a (l, 1')- We define (see p. 67 of Ref. [1]): 

A(p,R,T)= J dre^ r [ig<(r,t = 0;R,T)} 

= J dre-^ r (^(R - T -,T)MK + |, T)> (27) 

where, by definition, 

/ 7^/i(P' R > T ) =fi ( R > T )- (28) 

We can see that /i(p, R, T) corresponds to the well-known Wigner distribution function. In 
the classical limit, it reduces to the distribution function giving the number of atoms with 
momentum p at point R and time T. The symmetry-breaking terms in (6) result in the 
finite value of anomalous Green's function, and thus it is natural to introduce an additional 
distribution function for the non-condensate atoms which will give us the anomalous non- 
condensate density m(l), namely 

/ 2 (p,R,T)= J dre- ip - r [^<(r,t = 0;R,T)] 

= J dre-^iMK-^T^uiR+^T)}. (29) 

One can easily verify that the pair correlation function m in (16) is given by 

/ 7^/2( P ,R,T) = m(R,T). (30) 

It is important to remember that, as defined, the distribution functions fi and describe 
the behaviour of the atoms. They should not be confused with the distribution function for 
the quasiparticle excitations, such as discussed in Ref. Q| . 
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To obtain kinetic equations, we follow the Kadanoff-Baym approach in |1],|2| and rewrite 
the equations for ^(1,1') and g< 2 (l,l') in the relative and center-of-mass coordinates 
(r, t; R, T). We could obtain equations for g(r, t; R, T) as in Ref. 0, but for the simple 
HFB self-energies of interest it is sufficient to consider t=0, i.e. set ty = tf = T (see Ch. 7 
of [|J). Using (|24D , this procedure gives 

i-^= + MJ» ff (r, R, T) + -V r - V R + iv7(r, R, T) • V r + V + (r, R, T) ■ V R 



<9T 



m 



+ l -W r ■ v;(r, R, T) + 1 V R ■ v+(r, R, T) 



^n(r,R,T) 
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m(R + ^ T)^< (r, R, T) - m*(R - ^ T)p<(r, R, T) 



(31) 



where we have introduced the abbreviations 



AU» ff (l, 10 = //(l) - //(!') - I7 e// (1) + r/ e// (io 



(32) 



and 



v±(r, R, T) = v s (R + £ T) ± v s (R - £ T). 



(33) 



Similarly, the equation of motion for (?^(1, 1') expressed in the (r, R; £, T) variables gives 



+ /x(R + |, T) + //(R - J, T) - C/ e// (R + J, T) - [/ e// (R - £ T) + - 



4m 



+ m" V " + Ws " (r ' R ' T) ' Vr + ^ Vs+(r ' R ' T) ' Vr + \ Vr ' V7(r ' R ' T) 



+-V R -v+(r,R,T) 



#i 2 ( r > R ; T ) = 9 



m(R + |, T)^<(r, R, T) + m(R - |, T)g<(r, R, T) 



Finally, the equation of motion fl2q) for the amplitude of the order parameter can be written 
in the (R, T) coordinates as 

d dO(R,T) , V R 1 2 



<9T 



<9T 



2m 2 



mvf (R, T) + jj — U ext (R) - U(R, T) 



-g [2n(R, T) + n c (R, T)] + iv 8 (R, T) • V R + -V R ■ v s (R, T) 



n c (R,T) 



gm(R,T)Jn c (R,T)+ V '(R,T). 



(35) 
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The coupled set of HFB equations given by (pip , (|34J) and (|35| ) are the main formal 
results of this paper. These are a straightforward generalization of the analogous equations 
in the normal Hartree-Fock approximation discussed by Kadanoff and Baym (see Eq. (7.7) 
on pg. 71 of Ref. [[[]]). These results are important since they allow us to go beyond the 
simple HFP approximation which has been the basis of recent work on the non-equilibrium 
properties of a trapped Bose-condensed gas (see, for example, 0,0). 



We next proceed to use (|3lD and (plf ) to derive self-consistent equations for /x(p,R, T) 
and / 2 (p, R ; T) for the case when the external perturbation varies slowly in space and time. 
In this case, we expect that physical quantities v s , /i, U, U ext , n(l), and m(l) all vary slowly 
as functions of the center-of-mass coordinates (R, T). Thus, in the lowest approximation, 
and using the fact that small values of r are most important, we can use 



v s (R±^,T)=v s (R,T)± 



r _ 



v,(R,T), (36) 



and hence 



v+(r,R,T) ~2v s (R,T), vj(r,R,T) ~ [r • V R ] v s (R, T), 

AU? ff (v, R, T) ~ -r ■ V R [C/ e// (R, T) - //(R, T)\. (37) 

If we rewrite ([H|) using (|27|) and (|29|), and Fourier transform it, we obtain after some algebra 

^ - V R [e p + v s • p] ■ V p + V p [e p + v s • p] • V R ] / x (p, R, T) 
= -ig [m(R, T)/ 2 (-p, R, T) - m*(R, T)/ 2 (p, R, T)] 

+ | [V R m(R, T) ■ Vp/ 2 (-p, R, T) + V R m*(R, T) ■ V p / 2 (p, R, T)] , (38) 

where the "normal" single-particle energy is defined by 

e~ p (R,T) = |1 + U eff (R,T) - /i(R,T). (39) 

We emphasize that e p is not, in general, the local HFB excitation energy. The corresponding 
kinetic equation for / 2 (p, R, T) is: 
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^ + z2e p (R, T) - V R [v s ■ p] ■ V p + V p [v s • p] • V R ] / 2 (p, R, T) 
-igm{K, T) [A(p, R, T) + / x (-p, R, T) + 1] 
+|v R m(R, T) ■ V p [A(p, R, T) - f^-p, R, T)) . 



(40) 



These are the coupled HFB collisionless kinetic equations in a frame moving with the velocity 
v s , for the case of slowly varying disturbances. 

The equation of motion ( j35|) for the amplitude of the order parameter is an exact equa- 
tion, and is not limited for a slowly varying disturbances. Equating the real and imaginary 
parts of (|35l), we obtain two hydrodynamic equations of motion for the condensate 

<9n c (R,T) 



dT 

d9(R, T) 1 



•Vr K(R, T)v s (R, T)] + 2gn c (R, T)Im [m(R, T)] 



+ -mv s 2 (R, T)-fM> = -/i(R, T) - t/(R, T) 



5T ' 2 \»AR.T) 
where the condensate chemical potential yu(R, T) is defined by 



(41) 
(42) 



/i(R,T) 



VWn c (r,T) 



2mJn c (R,T) 



^(R) + t/(R, T) + g [2n(R, T) + n c (R, T)] + [m(R, T)] 



(43) 



Taking the gradient of (^2|), we obtain the generalized Landau equation 



<9v s (R,T) 1 9/T ^ mx 
V +2 VrV ' (R ' T) 



-V 



R 



fj,(R,T)+rj'(R,T) 



1 



(44) 



/n c (R,T)_ 

In ZGN 0, the contribution of the external symmetry-breaking field 7/ was left implicit. In 
the often used Thomas-Fermi (TF) approximation, the kinetic energy of the condensate is 
omitted and the HFB approximation for the condensate chemical potential then simplifies 
to 



/i TF (R, T) = ^(R) + 2gn(R, T) + gn c (R, T) + gRe [m(R, T)] . 



(45) 
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Next we consider Eq.(p8|) in the high temperature limit. In that case, n c is small and 
therefore we can neglect the terms proportional to gn c f2 in (|38| ) as small. Also, we know 
that m must be at least of order g |12[], and therefore the gift contribution to the self-energy 
is 0(g 2 ). Thus, we can neglect the right side of fl38|), leaving 



d 

~ V R [e p + v s • p] • V p + V p [e p + v s • p] • V R 



A(p,R,T) = 0. 



(46) 



This is precisely the expected collisionless Boltzmann equation for f\, valid in the HFP 
approximation. This approximation is only valid at finite temperatures, close to Tbec-, m 
which case e p (R, T) as defined in (|39|) is the correct excitation energy. In this limit, the 
kinetic equation fl46|) becomes equivalent to that derived (using a different formalism) by 
Kirkpatrick and Dorfman (KD) for a uniform gas. It is the local rest frame equivalent 
of the kinetic equation used by ZGN M if one ignores the collision terms. 



IV. STATIC HFB EQUILIBRIUM SOLUTIONS IN THE SEMI-CLASSICAL 

APPROXIMATION 



For a uniform system, the matrix Green's functions defined in (4) depend only on the 
relative coordinates, i.e. ^ Q/ g(l, 1') = g a /?(l — 1'). One can then solve (7)-(10) by Fourier 
transformation to obtain expressions for gn(p, oj) and g 12 (p, cu). Using these, we obtain the 
single-particle spectral density in the form M 



A u (p, w)= — (p, uj + iO H 



-21m 



uj - E p + i0+ u + E p + iO- 



2tt 



u 2 p 5(uj - Ep) - Vp5(u + E p 



Al 2 (p, oj)= —2Imgi 2 (p, uj + i0 + ) = 21m 



u p v p 



u p v p 



uj- E p + iO+ uj + E p + i0+ 



-2-KUpV* [5{uj - E p ) - 5{uj + E p )} . 



(47) 



The HFB excitation energy E p is given by 
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with 



and 



El 



V 

2m 



2gn — \i 



- (gmf 



e P=2m +29n ~ ^ 



u: 



\_E p 



UpV p 



gm 
2K, 



(48) 



(49) 



(50) 



The results in (f47|) have the same structure as in the simpler Bogoliubov approximation (see 
Ch.14 of Ref H). From (|17|), with $(1) = const., it follows that chemical potential in HFB 
approximation is given by MM 



fi = g {n + n + m) . 



(51) 



If we use this result in (EBI), it reduces to 



E 



p=0 



g 2 W 



~ 1 2 

m 



g [n c + m 



~ l2 



-Ag 2 mn c 



(52) 



Therefore, the long-wavelength HFB excitation spectrum has a finite energy gap. 

The density of non-condensate atoms in a uniform system can be found using the single- 
particle spectral density given in (]47| ) 

dp f duj 



n 



(2vr) 3 J 2tt 

dp 
(2^)3 

dp 
(27p 



N (lu)A u (p,uj) 
v 2 p +(u 2 p + vl)N (Ep 
^ [2N (E P ) + 1] - \ 



(53) 



where N (E P ) is the Bose distribution function. Similarly, the anomalous density is given 
by 
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m 



dp 



(2tt) 3 ^ p 
dp gm 



v p u* P [2N (E P ) + 1] 

[2N (E p ) + 1] . (54) 



(2vr) 3 2E p 

We want to find an approximate solution for f\ and f'2 for a trapped Bose gas which 
will be the analogue of the uniform gas results in (|53| ) and (|54"D. We start from the coupled 
static HFB equations, as derived in Ref. ||: 

Cui(R) - gm (R)vi(R) = ^Wj(R) 

Cvi(R) - gm (R)ui(R) = -E iVi (R), (55) 

where 

V 2 

£ = -— + U ext (R)-n + 2gn(R). (56) 



We can solve (pq) using the well-known semi-classical approximation. We assume that 
the normal and anomalous density are smooth functions of R on the scale of length of 



a,Ho = Y^^5 which defines the size of the condensate in a harmonic potential with a trap 
frequency u)o (we reinsert h in this discussion for physical clarity). Hence the self-consistent 
Hartree-Fock mean field varies slowly (as a function of R) on the length scale of order ano- 
Therefore, we can assume that Ui(R) and f«(R) have a form of a plane waves with a slowly 
varying amplitude in that region ||13|| , i.e. 



Ui(R) = u(p,R)e ilf ^ 

^(R)=^(p,R)e^ R ). (57) 

We introduce the momentum of elementary excitations by 

p = Wr^(R), (58) 
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which satisfy the condition p 3> h/ano- This condition, when expressed in terms of wave- 
length reduces to the small wavelength limit (A <C duo), or equivalently, to the semi-classical 
approximation limit expressed as kT ^> Hujq, where hu>o gives the harmonic well energy level 
spacing. If this condition is satisfied, we can neglect the spatial derivatives of u and v, as 
well as the second spatial derivative of ip. This is consistent with a so-called quasi- classical 
condition, which requires that a spatial change in wavelength of the particle must satisfy 
the condition d\/dx 1 [TJj]. The assumed form given by ( |57"D is only valid in the regions 



of space where this condition is satisfied. To treat the condensate in the corresponding 
approximation, we use the Thomas- Fermi approximation which is valid in the large N limit. 
The only region where the TF approximation for the order parameter is inadequate is close 
to the classical turning points at the condensate boundary [0,0 , which is consistent with 
inapplicability of the semi-classical approximation near these points. 

Putting all this together, we can easily solve the coupled equations (|55| ) for u(p, R) and 
u(p,R), 

gm R) 



2 / t> \ 6 p( r ) + „j/ ns e p [R)-E p [R) 



2E P (R) 



2E P (R) 



u(p, RK( P ,R) 



2E P (R)' 



(59) 



The local HFB quasiparticle energy E P (H) is given by 



E p 2 (R) = eJ(R) - (^mo(R)) 2 
We use the expressions for n and m given in terms of u and v 



- fio + U ext (R) + 2gn (R) 

2m 



1 2 



(gm (R)) 2 . (60) 



n(R) = ([l I 2 + I I 2 ] N (Ei)+ | Vi (R 

m(R) = -Du<(RK(R) [2iV (^) + 1] , 



(61) 



where N (E) is the Bose distribution for the quasiparticle excitations. The sum over the 
quantum states is replaced by the integral / dp/(27r) 3 and the semi-classical approximation 
for the diagonal distribution function / 10 (p, R) is given by 
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A„(p, R) = v'(p, R) + [« 2 (p, R) + v 2 (p, R)J N (E P (R)) 

while for the off-diagonal distribution function /2o(p, R) we obtain 

/2o(P ' R) = ~2E^j [2No{Epm + 1] • 



(62) 



(63) 



In summary, in the semi-classical approximation, the local static equilibrium normal 
density is given by 

dp 



no(R) = / 



(27T) 



^jTO(R) + l]-i 



(64) 



while the local static equilibrium anomalous density by 



m (R) = -prrio (R) / 



dp 

(2^ 



2jV (E p (R)) + 1 
2E P (R) 



(65) 



As expected, these semi-classical approximation results are the natural generalizations of 
the results obtained for a uniform gas given by fl53|) and fl54]) . The same kind of semi- 
classical results have also been obtained in Ref. [13[] for a trapped Bose gas using the Popov 
approximation (which corresponds to the HFB with m=0). Since the local quasiparticle 
energy E P (R) given by (|60|) depends on the normal and anomalous densities, the quantities 
in (|60|), (|64]) and ( |65]) must be solved self-consistently, as in Ref. fl~3[ . 



We can now show that these semi-classical static HFB results for fio and f'20 satisfy our 
collisionless static HFB kinetic equations. First of all, we note that in the static limit, fllCf ) 
reduces to 



2e p (R)/ 20 (p, R) = -gm {R) [/ 10 (p, R) + / 10 (-p, R) + 1] 

-i 9 - V R m(R, T) ■ Vp [/ 10 (p, R) - /io(-p, R)] ■ 



(66) 



In static equilibrium without any mass current, we have /io(p, R) = Ao( — P?R)> an d thus 
(|66|) can be further simplified to 



/20(P,R) 



2e p (R) 



[2/ 10 (p,R) + l]. 



(67) 



We note that the HFB approximation for e p (R) in (p9|) is given by 




(68) 



Similarly, in static thermal equilibrium, Eq.(^) for the diagonal distribution function 
/i(p, R, T) reduces to 



[-V R e p (R) • V p + V p e p (R) • V R ] /i (p, R) = sV R m (R) • V p / 20 (p, R). (69) 



In the last step, we again have used that /2o( — P, R) — /2o(P) R) an d that m is real. 

The results discussed above for /io(p, R) and /2o(p, R) describe local thermal equilibrium 
distributions functions induced by collisions, which are not included in our HFB collisionless 
equations. However, these local equilibrium functions do satisfy the static HFB kinetic 
equations. Inserting ( |B"2"D into (KBI), we obtain Q5"B|). This shows that the static equilibrium 
kinetic equation (|6~7D for the off-diagonal distribution function has a solution consistent with 
the semi-classical approximation. Substituting the local equilibrium distribution functions 
/io and /20 given by (|62f) and (^3|) , a lengthy but straightforward calculation shows that 
they satisfy the static HFB kinetic equation fl69"l). 

We recall that the semi-classical local equilibrium expressions for fi and / 2 are only valid 
under the condition ksT Tiujq and E P (R) 3> hu . This means that these approximate 
forms cease to be valid at very low temperatures (for further discussion, see Ref. [pT3| ). 



In this section, following the approach used in Ref. ||16|| , we prove that our coupled 



center-of-mass oscillation of the condensate and the non-condensate corresponds to 



V. KOHN MODE 



equations exhibit the rigid in-phase oscillation (or Kohn mode) ||. This rigid in-phase 
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n c (R,T)=n c0 (R-rj(T)) 
n(R, T) = n (R — r){T)) 

m(R,T) = m (R — rf(T)). (70) 

Here the center-of-mass displacement tj(T) (with v s = v n = r/) is independent of position 
and satisfies the harmonic oscillator equation of motion 

where uo a is the harmonic well trap frequency in the a th direction. 

We recall equation of motion for the superfluid velocity given by (ffij). For the case of a 
rigid in-phase oscillation described by (|70|), this equation can be written in the form 
d 2 r] 



m- 



-V R [/i (R - rt(T)) + ^(R) - U ext (K - r/(T))] . (72) 



Q T 2 

Since fiQ is position independent, it follows that the first term on the right hand side makes 
no contribution and we are left with 

= [f4rf(R) _ Uea *( R ' ^ ■ (73) 

For a harmonic trap potential described by 

U ext (R) = l -m{u 2 x x 2 + u 2 y y 2 + cu 2 z z 2 ), (74) 

(ff3|) reduces to (|7T|), as claimed. 

We next argue that the rigid in-phase oscillation described by ( |7D| ) corresponds to the 
following distribution functions: 

/ 1 (p,R,r) = / 10 (p,R-7 7 (T)) 

/ 2 (p,R,T) = / 20 (p,R-r7(T)), (75) 



where /i (p,R) and / 2 o(p?R) satisfy the static equilibrium kinetic equations (^) and ( |67D 
and the static equilibrium densities satisfy (^). In this case, (|38| ) for /i(p,R, T) reduces to 
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^ - V R [e p (R, T) + v s ■ p] • V p + V p [e p (R, T) +v s ■ p] ■ Vr 
= -ig [m(R, T)/ 20 (-p, R — 77) — m*(R, T)/ 20 (p, R - 77)] + 



/io(p,R- 77) 



+| [V R m(R, T) ■ V p / 2 (-p, R - tj) + V R m(R, T) ■ V p / 2 (p, R — 77)] . (76) 

From the definition of /jo(p, R), we see that f i0 (—p,H — r](T)) = /io(p, R — r](T)). We also 
note that for the in- phase oscillation under consideration, m(R, T) = mo(R — rj(T)), where 
mo is real. Therefore the first term on the right side of (|76|) vanishes, leaving us with 



^ - V R [e p (R, T) + v s • p] • V p + V p [ep(R, T) + v s ■ p] ■ Vr 
-- ^V R m (R - rj) ■ V p / 2 o(p, R - rf). 



/io(p,R- rj) 



(77) 



Using (|39D, ( ^3] ) and ([T0|), it is straightforward to see that e p (R, T) = e p (R — r]). Using this, 
([T?|) can be rewritten as 

" d 



dT 



- V R e p (R - r/) • V p + V p [e p (R - t]) + rj ■ p] • V R 



/io(p,R- rj) 



^V R m (R - 77) • V p / 20 (p, R - 77). 



(78) 



If we introduce new variable R'(T) = R — 77 (T), and note that 

d 

V r = Vr ^ = -77-Vr, 
then (|78|) can be rewritten as 



(79) 



[-77 • Vr - Vr6 p (R') ■ V p + V p e p (R) ■ Vr + 77 • Vr] /i (p, R') 

= 5 V R ™ (R').V p / 2 „(p,R'). (80) 



Since the first and the last term on the left side of (|80D cancel each other, we are left with 
an equation which is precisely the same as the static HFB equation (|69"D for the diagonal 
distribution function j\. 

Similarly, for the in-phase mode, Eq.(pE0|) for the off-diagonal distribution function f 2 
reduces to 
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^ + 2ie p (R ) T)+T 7 .V R 



/ 20 (p, R - rj) = -igm {H - 77) [2/ 10 (p, R — 17) + 1] . (81) 



Using (|79|), ( ]STD becomes 



/»(P, RO = "f^y [ 2 Ao(P, R') + 1] • (82) 

Again, this result is seen to be equivalent to the static HFB equation result given by ( |66[ ) 
and f3). 

In summary, we have explicitly verified that the ansatz given in ([TOD satisfies our cou- 
pled HFB kinetic equations for the two distribution function fi and fi- The oscillating 
center-of-mass displacement f}(T) satisfies the SHO equation of motion in ([H]). Thus, we 
have verified that our HFB coupled kinetic equations exhibit a solution corresponding to a 
rigid SHO oscillation of the static equilibrium density profiles (n c , n and m) in the direction 
of r) a , with the trap frequency u a . This expected solution is an important check on the 
correctness of our HFB equations of motion. We might note that we did not make use of 
the approximate semi-classical local equilibrium forms for f\ and / 2 given in Section IV to 
prove that the Kohn mode solution exists. 



VI. CONCLUSIONS 

To summarize, in this paper, we have used the HFB self-energy approximation to derive 
the equation of motion for the condensate order parameter (given by (17)) and the kinetic 
equations for the diagonal and the off-diagonal distribution functions for the non-condensate 
atoms (given by fl5S|) and (fiDD). We have only solved these coupled equations to exhibit the 
rigid in-phase Kohn mode (Section V). However, we emphasize that, more generally, these 
equations will lead to collective modes which satisfy various conservation laws ||, even 
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though they were generated from the HFB single-particle self-energies (which gives a single- 
particle spectrum with an energy gap in a uniform Bose gas). The essential physics is 



discussed in Sections 6 and 7 of Griffin's Varenna lectures |I2|. Another way of saying this 



is that for a uniform Bose gas, the coupled equations of motion in this paper will lead to 



the same "conserving" density response function discussed by Cheung and Griffin [|17| . 

Clearly the next step is to extend our present analysis by including the appropriate 
second-order self-energy contributions and hence, incorporate the effect of collisions into our 
kinetic equations. This generalization has been carried out for a uniform Bose-condensed 
gas by Kane and Kadanoff (KK) 0, work which has been recently extended to trapped Bose 
gases ||18||. KK concentrated on using their kinetic equations to derive the two-fluid hydro- 



dynamics equations. Kirkpatrick and Dorfman [|TTJ also derived a kinetic equation for the 
distribution function describing excitations in a uniform Bose-condensed gas, using a differ- 
ent formalism than KK. We believe that the KK approach based on non-equilibrium Green's 
functions gives the most systematic approach to deriving generalized quantum Boltzamnn 
equations for trapped Bose gases. 

ACKNOWLEDGMENTS 

We thank Tetsuro Nikuni and Eugene Zaremba for useful comments on the manuscript. 
This work was supported by a research grant from NSERC. M.I.T. would also like to ac- 
knowledge a graduate scholarship from NSERC. 



23 



REFERENCES 

[1] L.P. Kadanoff and G. Baym, Quantum Statistical Mechanics (Addison- Wesley, 1989). 

[2] J.W. Kane and L. P. Kadanoff, Journ. Math. Phys., 6, 1902 (1965). 

[3] H. Shi and A. Griffin, Phys. Rep. 304, 1 (1998). 

[4] A. Griffin, Phys. Rev. B 53, 9341 (1996) and references therein. 

[5] G.Baym and L.P. Kadanoff, Phys. Rev 124, 287 (1961). 

E. Zaremba, A. Griffin and T. Nikuni, Phys. Rev. A 57, 4695 (1998). 



[7] M.J. Bijlsma and H.T.C. Stoof, |cond-mat/9807051 . 



A. L. Fetter and J.D. Walecka, Quantum Theory of Many Particle Systems (McGraw-Hill, 
New York, 1971). 

[9] P.C.Hohenberg and PC. Martin, Ann. Phys. (N.Y.) 34, 291 (1965). 

[10] N.N. Bogoliubov, Lectures on Quantum Statistics, Volume II (Macdonald Technical and 
Scientific, London, 1971). 

[11] T.R. Kirkpatrick and J.R. Dorfman, Journ. Low Temp. Phys. 58, 308 (1985); 58, 399 
(1985). 

[12] A. Griffin, in Bose-Einstein Condensation in Atomic Gases, ed. by M. Inguscio, S. 
Stringari and C. Wieman (Italian Physical Society, Bologna, 1999). 

[13] S.Giorgini, L.P. Pitaevskii and S. Stringari, Journ. Low Temp. Phys. 109, 309 (1997). 

[14] L.D. Landau and E.M. Lifshitz, Quantum Mechanics, Ch. VII (Pergamon Press, Oxford; 
New York, 1977). 



24 



[15] F.Dalfovo, S. Giorgini, L.P. Pitaevskii and S.Stringari, Rev. Mod. Phys, in press; |cond- 



mat/ 9806038 . 



[16] E. Zaremba, T. Nikuni and A. Griffin, to be published. 

[17] T.H. Cheung and A. Griffin, Phys. Rev A 4, 237 (1971). 

[18] Milena Imamovic-Tomasovic, M.Sc. Report, University of Toronto (1998) 



25 



